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ABSTRACT 

Stars in open clusters are powerful probes of the intervening Galactic magnetic field, via back¬ 
ground starlight polarimetry, because they provide constraints on the magnetic field distances. We 
use 2MASS photometric data for a sample of 31 clusters in the outer Galaxy, for which near-IR po- 
larimetric data were obtained, to determine the cluster distances, ages, and reddenings via fitting 
theoretical isochrones to cluster color-magnitude diagrams. The fitting approach uses an objective 
minimization technique to derive the cluster properties and their uncertainties. We found the ages, 
distances, and reddenings for 24 of the clusters, and the distances and reddenings for six additional 
clusters that were either sparse or faint in the near-IR. The derived ranges of log(age), distance, and 
E{B — V) were 7.25-9.63, ^670-6160 pc, and 0.02-1.46 mag, respectively. The distance uncertainties 
ranged from ^8 to 20%. The derived parameters were compared to previous studies, and most cluster 
parameters agree within our uncertainties. To test the accuracy of the fitting technique, synthetic 
clusters with 50, 100, or 200 cluster members and a wide range of ages were fit. These tests recovered 
the input parameters within their uncertainties for more than 90% of the individual synthetic clus¬ 
ter parameters. These results indicate that the fitting technique likely provides reliable estimates of 
cluster properties. The distances derived will be used in an upcoming study of the Galactic magnetic 
field in the outer Galaxy. 

Subject headings: open clusters and associations: general—methods: numerical—magnetic fields: 
polarization—surveys—ISM: magnetic fields 


1. INTRODUCTION 


Open star clusters in the Galactic plane are useful tools 
to probe Galactic p roperties. Because cluster members 
are generally coeval (lTrumDleiilI925l : lFriell995ll . they are 
roughly the same age and located at the same distance. 
These properties make them excellent potential probes of 
the Galactic m agnetic fie l d via background starlight po - 
larimetry (e.g., maiHIT^ iHiTu^ri^ ISerkowskilh 96.511 . 
The measured linear polarization orientation of back¬ 
ground starlight traces the direction of the plane-of-sky 
component of the magnetic field located between the ob¬ 
server and the stars. The distance estimates provided by 
clusters can be used as upper limits to the proximity of 
the magnetic field along the line of sight. 

The overarching goal of the study is to explore the 
large-scale Galactic magnetic field in the outer Galaxy, 
especially how its behavior changes with distance and 
in the presence or absence of spiral arms. Using near- 
infrared (NIR) linear polarimetric observations of a sam¬ 
ple of star clusters, we probed the magnetic field toward 
the outer Perseus arm and its inter-arm regions. To 
find magnetic field properties in the interstellar medium 
(ISM) along these directions, however, we must establish 
cluster distances. In this work, we compare background- 


subtracted cluster color-magnitude dia grams (GMDs), 
creat ed using NIR 2MASS photometry (iSkrutskie et al.l 
l2006fl . to theoretical isochrones using a minimization 
approach to determine distances to the clusters, as well as 
their ages and reddenings. The cluster distances derived 
through this approach will be the basis for an upcoming 
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study of the Galactic magnetic field in the outer Galaxy 
(Hoq et al. 2016, in preparation). 

The cluster properties, including distances, of this sam¬ 
ple have been determined in previous studies. However, 
these studies used different datasets and different fitting 
techniques, and so have non-uniform systematics and un¬ 
certainties. In addition, uncertainties of the cluster prop¬ 
erties were often not reported. Therefore, we developed 
an approach to fit isochrones to establish cluster proper¬ 
ties and their uncertainties and applied it to a uniform 
dataset. 

Several recent studies, using optical and NIR photo¬ 
metric data sets, have developed objective methods to 
fit isochrones to cluster GMDs with the goal of estab- 
lish ing cluster param e ters, in cluding, but not limite d 
to. iNavlor fc Jeffriesl (1200^ . Ivon Hinnel et al.l (l2006fl. 
i Hernandez fc Valls-Gabaiid (1200811. iMonteiro et al ' 


( 201C b~ rMacieiewsk i fc Niedzielskil (12007), Alves et al 


(IM^. IDiasetalJ (1201211 ICurtis et al.l (12'oTl' 
r anes et al.l ( 20131 1. iLee et alJ (|2ni3[l . and iPerren et alJ 
( 2015fl . While all of these studies examined multiple 


clusters, their samples had little overlap with our sample 
of clusters (Section [5]). Therefore, we developed an 
objective method to find the properties of the present 
sample of clusters. 

The method presented here shares similar features to 
many of the studies listed above, such as searching a 
parameter space consisting of age, distance, and redden¬ 
ing by comparing stellar isochrones to cluster CMDs, 
but differs in other respects. For example, the present 
study does not employ a Bayesia n approach to search 
through the parameter space le.g.. IMonteiro et'mi2nifll : 


^ We also briefly explored the effects of varying metallicity as a 
parameter. 
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lAlves et 311120121 : iJanes et'al]l2013t iPerren et alj|201^ . 
but searches the complete parameter space. In addition, 
this technique allows the derivation of the cluster proper¬ 
ties without the need to identify which stars in the cluster 
field of view are cluster members and which belong to the 
field. This was essential, as establishing cluster member¬ 
ship selection based solely on photometric data can be 
challenging. 

The sample consists of 31 clusters in the outer Galaxy 
for which NIR polarizations were obtained. Seven clus¬ 
ters failed the initial cluster CMD-isochrone fitting pro¬ 
cedure because they were either too faint or too sparse. 
Of these seven, the distances and reddenings were suc¬ 
cessfully derived for six clusters by fixing their ages. For 
the 30 clusters that were fit, we found a median log(age) 
of ~9.2, a median distance of 2900 pc, and a median 
E{B — V) of ^0.5 mag. 

The structure of the paper is as follows. Section 2 
describes the cluster sample and the 2MASS NIR pho¬ 
tometric data. Section 3 describes the steps of the anal¬ 
ysis, the results returned by the fitting procedure, and 
the tests of the fitting technique. Section 4 discusses the 
distribution of the derived cluster properties and their 
comparisons to previous published results. Our findings, 
summarized in Section 5, establish distances with 8-20% 
uncertainties. 

2. CLUSTER SAMPLE SELECTION AND OBSERVATIONS 

To study the large-scale structure of the magnetic field, 
open clusters spanning a wide range of distance estimates 
and Galactic longitudes in the outer Galaxy were cho¬ 
sen for NI R polarimetric observ ations with the Mimir 
instrument (|Glemens et al.l[2007f) . The sample of 31 cho¬ 
sen clusters spans Galactic longitudes I ^119-232° and 
latitudes b —12 to -1-32°. Table [T] lists the 2MASS 
coordinates and key properties of each cluster. The clus¬ 
ters were originally optically-discovered and span only a 
modest range of reddening values. Some of the clusters, 
such as M 67, also lie well off the Galactic plane. 

Archival 2MASS J, H, and Ks photometric data were 
chosen to compare to theoretical isochrones for two main 
reasons: 1) 2MASS provides a consistent photometric 
dataset for the entire cluster sample, and 2) the all-sky 
coverage of 2MASS allows for comprehensive measure¬ 
ments of background stellar contamination for eve ry clus¬ 
ter (e.g.. lBonatto fc Bicall2007t lAlves et al.ll2ni^ . Using 
additional datasets, such as the Sloan Digital Sky Sur¬ 
vey (SDSS), might reduce the uncertainties of some of 
the derived parameters for some of the clusters, but such 
data were not available for all clusters in the sample. So, 
using these data for only some of the clusters would have 
introduced bias. 

We obtained 2MASS data from the IPAC Infrared Sci¬ 
ence Archive (IRSA) consisting of 50 arcmin diameter 
fields centered on each cluster. This size was chosen 
to sufficiently characterize the level of background stars 
even for clusters spanning large angular radii. (J — K) 
versus H CMDs were created from these data and the 
CMDs were used in fitting theoretical isochrones. The 
radius of the selected area was the same for each cluster, 
regardless of cluster size. 

3. FITTING THEORETICAL ISOCHRONES TO CLUSTER 

CMDS 


To find cluster properties, along with corresponding 
uncertainties, we developed a multi-step process to fit 
theoretical isochrones to cluster GMDs. These steps are 
summarized here and discussed in greater detail below. 
We used the PAdova and TRieste Stellar Evolution Gode 
PARSEC) isochron es from the Padova set of models 
Bressan et al.ll20ll) . Because the isochrones used are 
not continuous analytic functions, Monte Carlo realiza¬ 
tions of the isochrones were created for the purpose of 
comparing to cluster CMDs. 

For each cluste r, we created two Hess diagrams (e.g., 
lAlves et al.ll2012[l of the CMD densities of 2MASS stars: 
one of the region containing the cluster and one of a 
region used to characterize the background surrounding 
the cluster. The background CMD was subtracted from 
the cluster CMD, and Monte Carlo realizations of the 
CMD overdensity were created. Monte Carlo realiza¬ 
tions of the theoretical isochrone s , follo wing the initial 
mass function (IMF) of iChabrieil (120051) . were also cre¬ 
ated at discrete points in a parameter space of log(age), 
distance modulus [(m-M)/f], and color excess [E(J — AT)]. 
At each point, an isochrone realization was generated at 
that age and shifted by the (m-M)// and E( J — AT) val¬ 
ues. Each cluster background-subtracted CMD realiza¬ 
tion was compared to the set of isochrone realizations 
over the entire parameter space. A statistic was com¬ 
puted at each point in parameter space to determine the 
goodness-of-fit between each isochrone realization and 
each cluster CMD. The parameters of age, distance, and 
reddening that yielded the minimum for each cluster 
were adopted as best representing the cluster properties. 

3.1. Background Subtraction 

The contamination of cluster CMDs by foreground and 
background stars is less severe in the outer Galaxy than 
in the inner Galaxy. However, for some clusters in the 
sample, there remains a significant level of contamina¬ 
tion, especially for those close to the Galactic disk. To 
find accurate cluster properties, it is necessary to re¬ 
move the effects of contaminating stars. However, de¬ 
termining cluster membership on a star-by-star basis is 
difficult when using only photometric information. Sev¬ 
eral studies have used photometric properties, sometimes 
in conjunction with other available cluster data, to se¬ 
lect c l uster members (e g..lCurri e et al.||2 0inl : iMaia et HI 
120101: iDias et all 120121 : lAlves et al.ll2012t) . For the pur¬ 
pose of finding the best-fit isochrone for each cluster 
CMD while minimizing background contamination, we 
converted each cluster CMD into a CMD imag€0 from 
which a background CMD image could be subtracted. 

® Because there were several steps involved in the fitting proce¬ 
dure, a distinct name was given to the output of each step. We 
summ arize these names here for clarity. A “CMD image” (Section 
13.1.211 is a Hess diagram, created by binning the color-magnitude 
space of a cluster into pixels and counting the number of stars that 
fall into each pixel. A “CMD Poisson draw image” is created from 
a “CMD image” by drawing a value from a Poisson distribution 
for each pixel, where the mean of the Poisson distribution is equal 
to the number of stars in the pixel. A “cluster overdensity CMD 
image” is the result of subtracting an outer region CMD Poisson 
draw image fro m a co rresponding inner region CMD Poisson draw 
image (Section 13.1.311 . A “cluster Monte Carlo realization” is cre¬ 
ated by popul ating a cluster overdensity CMD image with discrete 
stars (Section 13.2t . An “isochrone realization” is a Monte Carlo 
realization consisting of individual stars, fo llowin g an IMF, that 
represents a given isochrone curve (Section 13. 3. 111 . A “synthetic 
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Fig. 1.— Radial density profile of the stellar counts of NGC 2099 
based on 2MASS photometric data, with a Gaussian fit overlaid as 
the solid black curve. The cluster radius of ~440 arcsec, defined as 
1.5x Gaussian width, is indicated by the dotted vertical line. The 
3.5 X Gaussian width location is outside the plot range. The mean 
background star count level of ^^5 counts arcmin”^ is indicated by 
the black dashed horizontal line. A King profile was also fit to the 
profile, shown as the dashed-dotted curve. 

These steps are detailed below. 

3.1.1. Background Stellar Densities and Cluster Extents 

Radial density profiles of the cluster star counts were 
created from the 2MASS photometric data. To establish 
the background stellar density and the cluster angular 
extent, the star counts were annularly binned, where the 
bins were 50 arcsec wide and centered on the stellar over¬ 
density in each region. Each cluster radial density profile 
was fit with a Gaus sian plus a uniform distribution (e.g., 
iMercer et al.ll200^ . 

For each cluster, the 50 arcmin diameter region was 
separated into three annular sub-regions: an inner re¬ 
gion, a buffer zone, and an outer region. We defined the 
cluster radius to be 1.5 times the Gaussian widtlQ (i.e., 
R-i. 5 Gw) and set the inner region radius to this Ri.sgw- 
The outer region was used to estimate the level of back¬ 
ground contamination of field stars in the inner region, 
and was defined as the area outside of 3.5 times the Gaus¬ 
sian width. The number of cluster stars, found by inte¬ 
grating the radial density profile, the background con¬ 
tamination level, and the Ri.sgw radius of each cluster 
are listed in Table [TJ 

The buffer zone was an annular gap, from 1.5-3.5 
Gaussian widths, located between the inner and outer 
regions. Stars located in this region were not used in 
the analysis. The use of the buffer zone was necessary 
to avoid regions where neither cluster members nor field 
stars dominate significantly, which would make it diffi¬ 
cult to separate the two groups. 

cluster” (Section l3.5ll is a cluster CMD realization generated from 
a theoretical isochrone at a specified age and metallicity, which is 
shifted by a given distance modulus and reddening, and whose re¬ 
alization stars were scattered in the CMD by typical photometric 
noise. 

The Gaussian width is equal to one Gaussian cr, but as we 
use the term a for the parameter uncertainties, the term Gaussian 
width is used here to reduce confusion. 


Figure [T] shows the radial density profile for cluster 
NGC 2099, with its Ri.sgw cluster radius indicated. For 
a rich, relatively nearby cluster, such as NGC 2099, the 
cluster angular extent is large. As can be seen from the 
figure, enlarging the cluster radius beyond Ri.sgw would 
likely introduce more background stars than cluster stars 
and increase the contamination of t he CMD. 

A King profile, following Eq. 14 of iKind (|1962h . which 
had the same number of terms as the Gaussian fit, was 
also fit to the density profile, as shown by the dashed- 
dotted curve in FigurelH While it appears to fit the den¬ 
sity profile well, the uncertainties returned (rcoj-e = 330 
± 70 arcsec) are larger than the uncertainties from the 
Gaussian fit (Ri.sgw = 448 ± 24 arcsec). In addition, 
it was not clear that a King profile was the appropriate 
function to fit to the sample, as some of the clusters are 
young [log(age)^7] and may not be dynamically relaxed. 
For these reasons, only Gaussian distributions were fit to 
the cluster radial density profiles. 

3.1.2. CMD Images 

For each cluster, separate {J — K) vs H CMDs were 
created for the inner and outer regions. For this process, 
and for the subsequent processes of creating cluster re¬ 
alizations and fitting to theoretical isochrones, we only 
used stars that were brighter than 15*^ mag in i7-band 
because the 2MASS photometric uncertainties increase 
significantly beyond this threshold. 

The individual stars in each GMD were binned into pix¬ 
els that were 0.1 mag wide in color and 0.25 mag wide in 
apparent magnitude to create stellar density GMD im¬ 
ages (Hess diagrams) that were 26x36 pixels. These pixel 
sizes were chosen to retain sensitivity to small variations 
in the star counts across the GMD, but also to be large 
enough to have significant star counts. The count value 
for each pixel was equal to the number of stars that fell 
into the pixel color and magnitude ranges. The outer 
region CMD counts were scaled by the ratio of the areas 
of the inner region to the outer region. Figure [5] (Top, 
A and B) shows the inner and outer region CMD images 
for NGC 2099. The inner CMD image (Panel A) shows a 
well-defined main sequence and a red clump, whereas the 
outer region image (Panel B) shows no cluster-like fea¬ 
tures. We conclude that the outer region contains very 
few cluster members. 

3.1.3. CMD Poisson Draw Images 

To create the cluster CMD image corrected for back¬ 
ground contamination, the simplest approach would be 
to directly subtract the area-scaled outer region CMD 
image from that of the inner region. However, doing so 
would yield incorrect uncertainties because stellar counts 
are governed by Poisson statistics in both regions. In¬ 
stead, a statistical approach was developed that pre¬ 
served the Poisson nature of the stars in all regions. 

A schematic of this approach continues in Figure [5] for 
NGC 2099. Panels A and B of Figure [5] show the inner 
and outer region CMD images. At every pixel in each 
of the CMD images, a value was drawn from a Pois¬ 
son distribution whose mean equaled the count in that 
pixeQ Panels C and D show the resulting CMD Poisson 
draw images of the inner and outer region CMD images, 

® For the outer region, a Poisson draw was created of the GMD 
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Fig. 2.— Illustration of the steps used to create an overdensity 
CMD image for NGC 2099. Top panels (A and B) show the inner 
and outer region CMD images, while the middle panels (C and 
D) show Poisson draws from the CMD images of A and B. The 
bottom left panel (E) shows the net overdensity of the star counts 
after the counts of D have been subtracted from those of C. Panel 
F shows the uncertainty of the overdensity image of Panel E. The 
size of each pixel in color-magnitude space is shown as a reference 
rectangle in the bottom right of each panel. The contour levels are 
drawn at [1, 5, 12, and 20] counts per pixel. 


respectively, for a single Poisson draw. These images 
closely resemble the CMD images of Panels A and B, 
but because they are draws from Poisson distributions, 
the counts of the pixels will fluctuate with each draw. 

The outer region CMD Poisson draw image was sub¬ 
tracted from the inner region CMD draw to create a clus¬ 
ter overdensity CMD image. Panel E of Figure [2] shows 
a cluster overdensity image for NGC 2099, which was 
produced by subtracting the Panel D outer region draw 
from the Panel C inner region draw. Panel F shows the 
Panel E image uncertainty, propagated from the Pois¬ 
son uncertainties of Panels C and D. A set of such clus¬ 
ter overdensity images for any one cluster will exhibit 
pixel count fluctuations governed by Poisson statistics 
that reflect the original inner and outer region properties. 
For sparse clusters, these overdensity images will exhibit 
more fluctuations because the lower cluster star counts 
will be more sensitive to background Poisson noise. 

NGC 2099 is a well-defined cluster, and therefore, is 
an easy case on which to perform the process illustrated 

image before the counts were scaled by the ratio of the areas be¬ 
tween the inner and outer regions. This Poisson draw was then 
scaled by the ratio of the areas. 


in Figure [5J We generated the same figure for King 1 
(not shown), which is a less populated cluster. The same 
features seen in Figure[2]for NGC 2099 are similarly seen 
for King 1, such as the prominence of the main sequence 
in the overdensity image, which indicates that the pro¬ 
cedure is robust. 

3.2. Creating Cluster Monte Carlo Realizations 

The method developed to match isochrones to clus¬ 
ter CMDs, described below, compares the colors and 
magnitndes of individual stars. Therefore, we converted 
each cluster overdensity CMD image into a cluster CMD 
Monte Carlo realization of individual points represent¬ 
ing stars. Each pixel of the overdensity CMD image was 
populated with a number of realization stars equal to the 
star count value in that pixel. The assigned colors and 
magnitudes of the realization stars in each pixel followed 
a uniform distribution within the bounds of the color and 
magnitude ranges of the pixel. Each realization star was 
assigned color and magnitude uncertainties based on the 
average 2MASS photometric uncertainties at that color 
and magnitude. These color and magnitude uncertain¬ 
ties were computed from the 2MASS photometric data 
of a randomly chosen 2.25 square degree field of view. 

To prevent outliers in the CMD Monte Carlo realiza¬ 
tion from adversely influencing the isochrone fit, stars 
far from the CMD main sequence and red clump were 
removed before comparing to isochrones. The Monte 
Carlo realization was temporarily converted to an image, 
similar to the process described in Section 13.11 by bin¬ 
ning the stars into color-magnitude pixels 0.1 mag wide 
in color and 0.25 mag wide in apparent magnitude. At 
each 0.25 mag step in magnitude, the number of stars 
per pixel along the color axis was fit with a Gaussian 
function. Individual realization stars that were located 
at colors further than three times the Gaussian width 
from the Gaussian peak were removed from the realiza¬ 
tion. No stars were eliminated when the Gaussian width 
was larger than 1 mag in color, as these magnitude bins 
were too sparsely populated to generate a significantly 
peaked overdensity. 

Figure [3] shows a GMD Monte Garlo realization for 
NGG 2099, which was fit by Gaussian functions of star 
counts along the color axis at each 0.25 mag step. Note 
that the actual binned pixel counts are not shown in the 
figure. The red clump and main sequence are well de¬ 
fined. Stars that are located further than three times the 
Gaussian width from the main sequence were flagged as 
outside the acceptable range and rejected from use (open 
black squares in the figure). The resulting “trimmed” 
cluster Monte Garlo realizations were used in the follow¬ 
ing steps. 

The number of stars populating the cluster Monte 
Garlo realizations fluctuated based on the Poisson draws 
of the inner and outer region GMD images. The exact 
color and magnitude locations of the realization stars 
within each color-magnitude pixel also fluctuated with 
each draw. Therefore, in the following procedure, we 
generated multiple cluster GMD realizations to capture 
the full ranges of these fluctuations. 

3.3. Isoehrone Fitting Procedure 

We developed a minimization appr oach to find the 
PARSEG isochrone (iBressan et al.l[201^ . at a given age 
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Fig. 3.— A Monte Carlo realization of an overdensity CMD im¬ 
age of NGC 2099. Triangles and squares represent realization stars. 
Centers of the Gaussian profiles of each magnitude step are indi¬ 
cated by red diamonds, and the [-Scr, H-Scr] width intervals of the 
steps are shown by the red horizontal error bars. Blue triangles 
represent realization stars that fall within three times the Gaus¬ 
sian widths. Black squares represent stars that fall outside this 
range, which were not included in the subsequent isochrone fitting. 


and metallicity and shifted by distance modulus and 
color excess, that best represented each cluster to find the 
cluster properties of age, distance, and reddening. Monte 
Carlo realizations of model stars were created from the 
isochrones. The cluster CMD Monte Carlo realizations, 
described above, were then compared to these isochrone 
realizations, as described below. 

3.3.1. Generating an Isochrone Realization 

Monte Carlo realizations of each isochrone were cre¬ 
ated following an IMF, where each isochrone was pop- 
ulated with indiv idual points representing stars (e.g., 

I Janes et al.l I2013I1 . These isochrone realizations were 
generated for every cluster, where the number of 
isochrone realization stars was set equal to the number 
of cluster CMD realizat ion stars. We ad opted a “kinked” 
IMF based on Eq. I in iChabried (I2005D : 


^{log (m)) = 


dn 


dlog (m) 


xexp 


= 0.093 


{log (to) - log (0.2))2 


0.605 
rUrmn <m< IMq 

= 0.041 , 7 ,,- 1 - 35 ± 0.3 ^rrimax > TO > IMq 


> ( 1 ) 


The mass range of each isochrone was uniformly di¬ 
vided into 50 bins, where the mass limits were set to the 
minimum (m,^i„) and maximum {m^ax) masses of the 
isochrone. These mass ranges were slightly different for 
every isochrone. Each mass bin was populated by the 
relative number density of stars per mass, based on Eq. 

[II 


Fig. 4. — An example of a Monte Garlo realization of an 
isochrone. The solid black line shows the PARSEG isochrone at an 
age of 1 Gyr and Z = 0.019. The blu e diamonds show 600 Monte 
Garlo realization stars that follow the lGhabrieil II2005I) IMF. 


Figured] shows a 1 Gyr, solar metallicity (J—K) versus 
H PARSEC isochrone with a Monte Carlo realization 
overlaid. The realization contains 600 stars, shown as 
open blue diam onds, which are distributed following the 
IChabrieil (I2005D IME. The isochrone realization reveals 
the expected concentration of stars at the theoretical red 
clump location and along the main sequence. 

3.3.2. Finding the Best-Fitting Isochrone 

To find the best-ht isochrone for each cluster, we tested 
isochrones across a 3D parameter grid of log(age), dis¬ 
tance modulus [(m-M)//] (uncorrected for extinction), 
and color excess [E(J — K)]. At each grid point, an 
isochrone was selected at that age and shifted in magni¬ 
tude by the distance modulus and shifted in color by the 
color excess. 

We tested whether having metallicity as another free 
parameter affected the results of the fitting procedure, 
and f ound that, in most cases, it did not (see Section 
I3.4.I|) . Based on these results, metallicity was not fit for 
the majority of the clusters, but was instead fixed for 
each cluster. If metallicity had been determined spectro¬ 
scopically for a cluster in a previous study, that metal- 
licity was adopted (see Table [T]). Otherwise, solar metal- 
licty was assumed. The exception was cluster NGC 1857, 
for which a metallicity of Z = 0.005 was used. In this 
case, solar metallicity isochrones would not fit both the 
main sequence and the small apparent red clump. We 
attempted to fit isochrones using several different metal- 
licities, and determined that a metallicity of Z = 0.005 
provided the best fit to the cluster CMD. 

The isochrone metallicity used for each cluster is listed 
in Tabled) While the PARSEC isochrones adopt a value 
of Z = 0.0152 for solar metallicity, to remain consistent 
with earlier studies, we adopted a metallicity of Z = 0.019 
(e.g.. lAnders fc Grevessell9^lGirardi et al.ll2000D as the 
solar metallicity. Therefore, in the context of the PAR- 
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SEC models, we have selected a solar metallicity that is 
somewhat metal-rich. 

The 3D parameter grid consisted of 20 steps in each di¬ 
rection of log(age), distance modulus, and color excess. 
We performed the isochrone fitting procedure, described 
below, twice: once for a coarsely stepped grid of param¬ 
eters and a second time for a more finely stepped grid. 
For the coarse grid, log(age) was stepped by 0.05, dis¬ 
tance modulus by 0.075 mag, and color excess by 0.025 
mag. The initial center of each cluster’s coarse grid in 
parameter space was determined by visually identifying 
the isochro ne that appeared to best match the cluster 
CMD fe.g.. lAlves et al.ll2012D . The finer grid, centered 
on the best-fitting point found from the coarse grid, con¬ 
sisted of log(age) steps of 0.025, distance modulus steps 
of 0.0375 mag, and color excess steps of 0.0125 mag. 

A cluster CMD Monte Carlo realization, generated fol¬ 
lowing the steps in Section [3.21 and an isochrone realiza¬ 
tion, following Section 13.3.11 with the same number of 
stars as the cluster realization, were compared at each 
parameter grid point. One of the key goals of the fitting 
procedure was to give more weight to stars with lower 
photometric uncertainties, i.e., brighter stars. Therefore, 
cluster CMD realization stars were ranked by bright¬ 
ness. A scaled distance was calculated between 

the brightest cluster star and each isochrone realization 
star in CMD space, defined as: 

CT(J-iC)^ CTHc 

( 2 ) 

where C refers to the cluster realization star and I refers 
to one isochrone realization star. {J — K)c and {J — K)i 
are the colors of the cluster star and the isochrone star, 
respectively. He and Hj are the H magnitudes of the 
cluster star and isochrone star, respectively, and 
and a He the uncertainties of the {J — K) color and H 

magnitude of the cluster star, respectively. Because the 
cluster stars came from Monte Carlo realizations, they 
do not have intrinsic uncertainties. So, the uncertainties 
assigned to them were based on the average 2MASS un¬ 
certainties at the same J, i7, and K magnitudes as the 
cluster stars. 

The isochrone star that yielded the lowest Dsj value 
was paired to the cluster star and removed from the 
pool. Another set of Dgj- values was computed using the 
next brightest cluster realization star and the remaining 
isochrone realization stars, and the isochrone star that 
yielded the lowest Dg^ value was paired to the cluster 
star. This process was continued until every cluster star 
was paired to one isochrone star. 

A representation of this process is shown in Figure 
[5l Select cluster realization-isochrone realization stellar 
pairs are highlighted. The brightest cluster realization 
star, at ^7 mag in iJ-band, is paired to the isochrone 
star that yielded the lowest Dsj value. Each cluster re¬ 
alization star was paired with an isochrone star, but only 
10 cluster stars, every 60th in the brightest ranked star 
list, are highlighted in Figure [SI Allowing the brightest 
clusters stars to be matched first ensured that these stars 
had somewhat greater influence on the fit. 

Once every cluster-isochrone stellar pair was created, 
the Ds, values of every pair were squared and summed to 



J-K (mag) 

Fig. 5. — Example of the process of matching the cluster CMD 
Monte Carlo realization stars to the isochrone realization stars to 
determine the best-fit isochrone for NGC 2099. The '^630 clus¬ 
ter realization stars are shown as black triangles and the ~630 
isochrone stars are shown as blue diamonds. Ten cluster stars (ev¬ 
ery 60th star, ranked by brightness) are shown connected to their 
paired isochrone stars with red lines, where the paired cluster and 
isochrone stars are displayed as red triangles and green diamonds, 
respectively. Note that this example represents a poor match be¬ 
tween the isochrone and cluster CMD, and was chosen to better 
display the pairing process. 

create a statistic to represent the goodness-of-fit be¬ 
tween the cluster Monte Carlo realization and isochrone 
Monte Carlo realization: 

= ( 3 ) 

p 

where p represents each cluster-isochrone stellar pair. 

For every cluster, 30 sets of Poisson draws were created 
from the inner and outer region CMD images, from which 
30 CMD overdensity images were created. One Monte 
Caro realization was created from each of the 30 overden¬ 
sity images, resulting in thirty cluster CMD Monte Carlo 
realizations for each cluster. Ten Monte Carlo isochrone 
realizations were created at each point in the 3D param¬ 
eter grid. Therefore, at each point in the parameter grid, 
to assess a particular isochrone’s goodness-of-fit, 300 
values were computed via Eqs. H] and |31 These 300 
values were median-filtered and a mean was computed 
to represent the goodness-of-ht at each grid point. This 
is a more conservative approach than simply selecting 
the smallest of the 300 values, which could easily be 
an outlier. 

The point in the 3D parameter grid that had the lowest 
mean x^ was chosen as the point whose isochrone param¬ 
eters of log(age), distance modulus, and color excess best 
represented the actual cluster properties. 

The uncertainties of the best-fit parameters were esti¬ 
mated by evaluating the mean x^ values of the grid. To 
estimate the uncertainty in log(age), the distance mod¬ 
ulus and color excess were held constant at their best-fit 
points, and the mean x^ values of the 20 points along the 
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Fig. 6.— NGC 2099 overdensity CMD image with the best-fit 
isochrone overlaid. Cluster overdensity contours at levels of [1, 
5, and 12] counts per pixel are shown in black, and the best-fit 
isochrone is shown as the solid blue line. A Monte Carlo realization 
of the best-fit isochrone is overlaid, with ^^630 stars shown as red 
diamonds. The best-fit parameters are listed, along with their 
corresponding uncertainties. The best-fit isochrone does not fall 
directly in the center of the red clump (see text). This isochrone 
was shifted by ±l(j in E(J — K) and drawn as the dashed green 
lines. The size of each pixel in the CMD image in color-magnitude 
space is shown as a reference rectangle in the bottom right of the 
panel. 


log(age) axis of the parameter grid were selected. The 
uncertainty in log(age) was found by calculating the 
weighted deviation of these 20 log(age) values: 


- XminYIX^) 


E(i/xD 


(4) 


where the Xi are the 20 log(age) values along the log(age) 
parameter axis and the xf represent the 20 correspond¬ 
ing x^ values. The best-fit value of each parameter is 
then Xmin ± where x is either log(age), distance 

modulus, or color excess. 

One cluster overdensity CMD image for NGC 2099 is 
shown in Figure IHl with the best-fit isochrone overlaid. 
This isochrone does not appear to fit the red clump per¬ 
fectly, but it does fall directly down the middle of the 
main sequence. This discrepancy may be due to the rel¬ 
atively small number of stars in the red clump compared 
to the large number along the main sequence, which dom¬ 
inates the fit. As can be seen by the ±lcr E{J — K) 
isochrones overlaid on the plot (dashed green lines), the 
— Icr limit isochrone that falls in the middle of the red 
clump does not fit the main sequence as well. This re¬ 
sult indicates that the isochrone selected does best rep¬ 
resent the overall shape of the cluster CMD at the given 
spectroscopic metallicity. We also tested whether using 
different metallicities would affect the fit, and those re¬ 
sults are described in Section [3.4. II 


3.4. Results of Isochrone Fitting 

Of the 31 clusters in the NIR polarization-based sam¬ 
ple, 24 were fit successfully. We determined whether a 


cluster was successfully fit by a visual comparison be¬ 
tween the best-fit isochrone and the cluster CMD. For the 
seven that were not fit successfully, the best-fit isochrones 
failed to overlap the cluster main sequence and/or the red 
clump of the CMDs. These seven clusters, Berkeley 14, 
Berkeley 18, Berkeley 70, Basel 11b, Berkeley 32, NGC 
2126, and Berkeley 39, were either faint in the NIR or 
sparse, thereby limiting the success of the fitting proce¬ 
dure. 

Figure [7] shows best-fit isochrones for the 24 
successfully-fit clusters. While the clusters are listed in 
order of Galactic longitude in the tables, they are ordered 
by age in Figure [7] to show the progression of the shape 
of the cluster GMDs and the distribution of stars as a 
function of age. As clusters age, the red clump becomes 
more populated. As they age further, the asymtotic gi¬ 
ant branch (AGB) becomes more populated as more stars 
evolve off the main sequence. This evolution of the loci 
of stars in cluster GMDs can be seen in Figure [71 The 
youngest cluster, NGC 869 has no observed red clump, 
whereas the intermediate age NGC 2099 contains a red 
clump. In the oldest cluster, Trumpler 5, the evolved 
stars are located along the AGB. 

Figure [5] is similar to Figure [71 but instead of plotting 
the best-fit isochrones on the cluster overdensity GMD 
images, the isochrones are plotted on GMDs of the indi¬ 
vidual 2MASS stars from the inner region of each clus¬ 
ter. Because all stars located in each inner region are 
included, the contamination from background stars is 
present. Nevertheless, the best-fit isochrones trace the 
loci of cluster stars. 

For the seven clusters that were not fit successfully 
for all three parameters, we repeated the procedure de¬ 
tailed above, but kept ages fixed, and only searched the 
2D parameter space of distance modulus and color ex¬ 
cess. The fixed age of each cluster was adopted from 
the most recent publishe d study to derive the properties 
([Kharchenko et ^l2013[) of these clusters. By only fit¬ 
ting for two parameters, we recovered the distance modu¬ 
lus and color excess for six of the seven clusters. Berkeley 
70 was still not properly fit, as it is too faint. Figured 
shows the best-fit isochrones overlaid on the six cluster 
overdensity GMD images and the 2MASS stars from the 
inner regions, similar to Figures |7| and [S] 

Table [2] lists the best-fit parameters of log(age), (m- 
M)//, and E( J — K), along with their corresponding un¬ 
certainties for the 30 clusters that were fit. The distances 
(in pc) and E(i? — V), derived from the (m-M)// and 
E(J — AT), are also listed. E(B — V) was calculated by 
equating Ay = RvE{B—V) and Ay = rE{J—K), where 
E{B -V) = rE{J-K)/Ry Ki E{J - K) /0.53. Ry and 
r were assumed to be 3.1 and ^5.9, respe ctively, appro¬ 
priat e for diffuse regions in the ISM (e.g., iWhittet et al.l 

[2?iol . 

The approximate numbers of stars in the cluster GMD 
Monte Carlo realizations are also listed in Table m The 
numbers of CMD stars were estimated by subtracting the 
total (area-scaled) star counts of each outer region CMD 
image from that of its inner region CMD image, which 
was essentially the total number of counts in the cluster 
CMD overdensity image. The uncertainties of the num¬ 
ber of stars were propagated from the star counts of the 
CMD images of the inner and outer regions. The num¬ 
ber of stars in each CMD Monte Carlo realization will 
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Fig. 7.— Similar to Figure [6l but for the 24 clusters for which all three parameters of log(age), distance modulus, and color excess were 
fit. Contour levels are shown at [1, 2, 5, and 7] stars counts per pixel. The clusters are plotted in order of increasing age to show the 
evolution of the shape of the isochrone as well as the distribution of stars in the CMD. 


fluctuate about this number based on the Poisson draw 
from each cluster CMD overdensity image. These Table 
[2] numbers are less than the number of cluster stars esti¬ 
mated by integrating over the cluster radial density pro¬ 
files (listed in Table [T]). This is because Table [2] only ac¬ 
counts for stars within the cluster inner regions (Ri.sgvf) 
and does not include stars fainter than 15*^ mag in H- 
band. 

3.4.1. Testing Isochrone Fits For Metallicity Dependence 

To determine whether adding metallicity as a free pa¬ 
rameter would result in improved isochrone fits to the 
cluster CMDs, we repeated a part of the fitting proce¬ 
dure for five different metallicities for five of the clusters. 
We fit the fine grid of parameters at metallicities of Z = 
0.005, 0.008, 0.012, 0.017, and 0.019. The five clusters 
selected were King 1, NGC 1245, King 7, NGC 2099, and 
NGC 2141. These were chosen because they had a vari¬ 
ety of ages, and had a combination of spectroscopically 
determined (NGC 1245, NGC 2099, NGC 2141, as listed 
in Table 1) and assumed (King 1, King 7) metallicities 


as noted in Table [TJ 

For each of the five clusters, the resulting five minimum 
values of the best-fit locations in the parameter grids, 
corresponding to the five metallicities, either remained 
unchanged or increased at low metallicities compared to 
the minimum value of the metallicity initially used in 
the isochrone fits, listed in Table [TJ For all five clusters, 
the three best-fit parameters found at each metallicity 
were within Ict of the best-fit parameters listed in Table 
2. Therefore, varying the metallicity did not change the 
best-fit parameters by more than la. Visual inspection 
of the best-fit isochrones overlaid on the cluster CMDs 
also showed no improvement, with the exception of NGC 
2099. At a metallicity of Z = 0.008, the red clump of the 
best-fitting isochrone was shifted blueward by about 0.03 
mag in (J — K) compared to the red clump location of 
the Z = 0.019 best-fit isochrone (shown in Figure |6|). 
While this shift places the isochrone red clump closer to 
the center of the cluster red clump, the amount of the 
shift is well within the fitting uncertainties of the color 
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Fig. 8.— Similar to Figure [T] but with the individual 2MASS (J — K) and H colors and magnitudes of all stars located in each cluster 
inner region plotted as black diamonds. The best-fit isochrones trace the cluster features for the 24 clusters. 


excess and is not statistically significant. Therefore, the 
best-ht parameters found with Z = 0.019 for NGC 2099 
are reported in Table [2] We conclude that the NIR colors 
are not ideal for determining metallicity, as optical colors 
are more sensitive to metallicity changes. 

3.5. Quality Testing using Synthetic Clusters 

To determine the accuracy of the fitting procedure, we 
created and fit synthetic clusters using the same proce¬ 
dure as described above. The synthetic clusters were 
generated from isochrone realizations at 12 ages, ranging 
from log(age) of 7.3 to 9.6, created using the procedure 
described in Section 13.3.11 The {J — K) colors of the 
stars were shifted by 0.25 mag, and the distance modu¬ 
lus (uncorrected for extinction) by 10 mag, to simulate 
clusters with reddenings of E(i3 — V) of 0.47 mag at 
distances of '^900 pc. The stars were then distributed 
about their color and magnitude values based on the av¬ 
erage 2MASS uncertainties corresponding to those colors 
and magnitudes. Clusters consisting of such “synthetic” 
stars at the 12 given ages were generated with each of 


50, 100, and 200 cluster members. 

Background held contamination was created using 
2MASS photometric data from a 2.25 sq. de gree h eld. 
Similar to the procedure described in Section 13.1.21 an 
outer region CMD image was created from this held. 
Each synthetic cluster was assigned a “cluster” radius 
of 160 arcsec, the average of the cluster radii listed in 
Table m to scale the counts of the outer region CMD 
image. 

The synthetic clusters consisting of individual stars 
were converted to CMD images, and then added to Pois¬ 
son draws of the outer region CMD image to simulate 
inner region CMD images (a cluster with held contami¬ 
nation). Poisson draws were created from the inner and 
outer region CMD images. The Poisson draws of the 
outer region images were subtracted from those of the 
inner region to create CMD overdensity images. Monte 
Carlo realizations of the CMD overdensity images were 
compared to Monte Carlo realizations of isochrones in 
3D parameter grids of log(age), distance modulus, and 
color excess. Thirty Poisson draws of the inner and outer 
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Fig. 9. — Similar to Figures[7]and|8]but for the six clusters where 
log(age) was fixed and only distance modulus and color excess were 
fit. The black contour levels at [1, 2, 5, and 7] counts per pixel show 
the cluster overdensity CMD images, and the black diamonds are 
individual 2MASS stars located in the cluster inner regions. 


regions were created, from which 30 overdensity images 
were created. One Monte Carlo realization was created 
from each overdensity image, resulting in 30 Monte Carlo 
realizations for each cluster. Ten isochrone realizations 
were created at each grid point. The parameters of the 
grid point that yielded the lowest mean value were 
adopted as most representative of the cluster properties. 

In total, 36 synthetic clusters were fit, consisting of 
three different numbers of members at 12 different ages, 
all with constant distance and reddening. 

Figure [10] shows the resulting best-fit age as a func¬ 
tion of the synthetic cluster input age, where the dashed 
lines represent offset lines of equality. Most of the fit 
ages (>88%) fall within their Icr uncertainties of the in¬ 
put ages. The largest deviations from the line of equal¬ 
ity are seen in the clusters with 50 stars. This result 
is reasonable, as these clusters have the least number of 
stars with which to define their main sequences. They 
are also the most sensitive to Poisson fluctuations and 
background contamination. The points deviate the most 
around log(ages) of 9.2 - 9.5 as the AGB is not well pop¬ 
ulated for clusters that have so few stars. 

Figure |TT] (A, top) shows the distance modulus re¬ 
turned from the isochrone fitting as a function of input 
age. The input distance modulus (10 mag in H-band) 
is marked by the dashed lines (with offsets). Similarly, 
for E(J — K), shown in Figure [TT] (B, bottom), the in¬ 
put E( J — K) of 0.25 mag is indicated by dashed lines. 
For both distance modulus and color excess, larger devi¬ 
ations are seen mostly at the youngest ages. We suspect 
these deviations are due to the absence of a red clump 
in the cluster CMDs, which leaves the fit somewhat un¬ 
constrained. More than 85% of the fit distance modulus 
values, and all the fit color excess values, fall within Icr 
of their input values. 

Based on these analyses of fitting to synthetic clusters, 
for which we know the true parameters a priori, we con¬ 
clude that the overall isochrone fitting procedure is re¬ 
liable within the uncertainties it returns, provided that 
the cluster membership is not overly sparse. Addition¬ 
ally, the prevalence of deviations that are less than Icr 



Fig. 10.— Input log(age) values of synthetic clusters are plotted 
against output best-fit ages for clusters with 50 (black diamonds), 
100 (blue triangles), and 200 (red squares) member stars. The 
black, blue, and red dashed lines show where input log(age) = fit 
log(age) for 50, 100, and 200-star clusters, respectively. The 100- 
star and 200-star cluster dashed lines have been displaced vertically 
by 0.4 and 0.8, respectively, for clarity. Most of the fit ages are 
equal to the input age within the fit uncertainties. Deviations are 
largest for the 50-star clusters, especially at the older ages. 
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Fig. 11.— The best-fit distance modulus (A, Top, uncorrected 
for reddening) and E(J — K) (B, Bottom) of the synthetic clus¬ 
ters are plotted against input log(age) of the synthetic clusters. 
The 50-star, 100-star, and 200-star clusters are shown as black di¬ 
amonds, blue triangles, and red squares, respectively. The 100-star 
and 200-star clusters were displaced vertically by 0.8 and 1.6 mag, 
respectively, for distance modulus, and by 0.2 and 0.4 mag, respec¬ 
tively for color excess. The dashed lines show the input parameters 
of 10 mag for (m-M)if and 0.25 mag for E(J — K). 


may indicate that the calculated uncertainties are some¬ 
what overestimated. We elect to retain them, though, as 
conservative uncertainty estimates. 

4. DISCUSSION 
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To better understand the properties of the cluster sam¬ 
ple, we searched for trends among the cluster parameters. 
Next, we compared the results of the present study to 
those found in previous studies to determine whether, 
and to what degree, our findings differ from previous 
ones. 

4.1. The Properties of the 30 Clusters fit to Isochrones 

The properties of the 30 clusters, determined via fitting 
theoretical isochrones, span wide ranges of age, distance, 
and reddening. The standard deviations of the distribu¬ 
tions of the best-fit log(age), distances, and E{B — V) are 
0.63, 1400 pc, and 0.3 mag, respectively, and the median 
values are ~9.2, 2900 pc, and 0.5 mag, respectively. The 
farthest cluster is Berkeley 60, at just over 6 kpc, and 
the nearest cluster is M 67 at ^670 pc. 

The subsample of six clusters for which only distance 
modulus and color excess were fit have relatively larger 
distance and reddening uncertainties than the subsample 
of 24 clusters for which all three parameters were fit. 
The mean uncertainties for the six clusters were ^^420 pc 
and 0.15 mag for distance and E(B-V), respectively. For 
the 24 clusters, the mean uncertainties were ~320 pc, 
0.14 mag, and 0.13 for distance, E(B-V), and log(age), 
respectively. 

The six clusters have somewhat larger parameter un¬ 
certainties due to their sparse and/or faint nature. As 
can be seen in Figure IHl Basel 11b and NGC 2126 have 
few cluster members in their CMDs, ^30 stars, which 
may not be enough to reliably constrain the fit. The 
four clusters Berkeley 14, Berkeley 18, Berkeley 39, and 
Berkeley 32 all have main sequence turnoffs near ^fi¬ 
ll.5 mag in i7-band, just above our H = 15 mag limit. 
This does not provide an adequate portion of the cluster 
main sequence to be fit. 

Based on the numbers of CMD stars of the clusters that 
could not be fully fit, and from the isochrone fits to the 
synthetic clusters, we expect that the fitting procedure is 
less robust for clusters with fewer than ^50 stars in their 
CMDs or whose main sequence turnoff is below H ^14 
mag. For these faint clusters, deeper photometric data 
would enable more successful fits. 

We note that NGC 2266, which was part of the 24 
successfully fit subsample, has fewer than 50 stars in its 
CMD. Some caution may be advised when adopting its 
parameters from this study. 

4.1.1. Comparison of Derived Cluster Properties 

The derived parameters of the 30 clusters were 
searched for correlations for potential biases in the fitting 
procedure. Figure[T3plots the best-fit log(age), (m-M)//, 
and E( J — AT) values against each other and against the 
number of CMD members. A linear fit was computed 
for each comparison, and the slope of each fit (m), along 
with its reduced value, are listed in each plot. 

No obvious trends exist between the number of CMD 
stars and the three fit parameters. There is a corre¬ 
lation between the distance modulus and color excess, 
which is reasonable given that a cluster is likely to be 
more extincted if it is farther away. Both distance mod¬ 
ulus and color excess appear to decrease as a function 
of increasing cluster age. The trend of distance modu¬ 
lus decreasing with age may be because younger clusters 



Fig. 12.— Best-fitting parameters of the 30 clusters whose CMDs 
were reliably fit, plotted against the other parameters to search 
for trends. (Left column.) Log (age), (m-M)^ , and E(J — K) 
plotted against number of CMD stars. (Middle column.) (m-M)^ 
and E(J — K) plotted against log(age). (Right column.) E(J — 
K) plotted against (m-M)ff. Linear fits to each set of values are 
overlaid as black solid lines, and the slope (m) and reduced 
value of each fit are reported. 

contain relatively brighter stars that can be seen at larger 
distances. The correlation between age and color excess 
may be due to the likelihood o f finding old er clusters far¬ 
ther from the Galactic plane (jFriell [l99^ . which would 
be along less extincted sight-lines. 

4.2. Comparison to Previous Studies 

We compared our derived cluster parameters to those 
found in recent publications to determine whether there 
were significant differences. Table [3] summarizes the rel¬ 
evant parameters found in the previous studies for the 
30 clusters. The fifth column of the Table lists whether 
our cluster log(ages), distances (in pc), and reddenings 
(E{B — V)) agree with those in the cited sources to within 
Icr, 2cr, 3cr, Or > 3(7 (a, b, c, d labels, respectively), 
based on our uncertainty estimates for each parameter. 
Log(age), distance in pc, and E{B — V) values were used 
for comparisons because these were the most cited prop¬ 
erties in the literature, especially as most of the studies 
were done in the optical. We find that for most of the pa¬ 
rameters of the majority of clusters (^80% of all the clus¬ 
ter parameters), the derived values agree to within 2cr. 
As no consensus of concordance was found in the litera¬ 
ture for comparisons of the agreement of derived cluster 
properties, we adopted 2 ct as the standard of agreement. 

Figure [T31 plots the distributions of the differences be¬ 
tween the parameters derived in this study to those of 
previous studies. No systematic offsets were found be¬ 
tween our values of log(age), distance, and reddening and 
those cited. Reddening estimates show the most agree¬ 
ment (> 95%) with previous values, while the derived 
ages and distances agree to within our uncertainties for 
>85% of the cited values. 

4.3. Galactic Locations and Cluster Properties 

Figure [14] plots the locations of the 30 fit clusters, with 
the Sun at the origin, using the derived distances listed 
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Fig. 13.— The distributions of differences between the cluster 
parameters derived in the present study and those derived in pre¬ 
vious studies. The blue striped distributions indicate the numbers 
of comparisons that differed by greater than 2(T. One point in the 
age distribution, four points in the distance distribution, and two 
in the K{B — V) distribution are not shown, as they are outside 
the plot ranges. 

in Table |2| The Perseus Spiral Arm (|Reid et al.l I2014H 
is denoted as the grey stripe. The clusters span a wide 
range of distances, residing in both arm and interarm lo¬ 
cations. The clusters located beyond the Perseus Arm 
have slightly larger reddenings on average than the clus¬ 
ters located in front of the arm. These distances will 
be combined with NIR polarimetry of the cluster stars 
to probe the nature of the magnetic field in the outer 
Galaxy. 


5. SUMMARY 

To study the large-scale structure of the Galactic mag¬ 
netic field in the outer Galaxy, we obtained NIR stel¬ 
lar polarimetric measurements of a sample of 31 open 
clusters spanning wide ranges of longitude, distance, and 
reddening. It was essential to determine the cluster dis¬ 
tances accurately to optimally use the polarimetric in¬ 
formation they provided. To do so, we developed a 
minimization technique to fit theoretical isochrones to 
the cluster GMDs. These fits returned the cluster prop¬ 
erties of distance, age, and reddening, along with their 
corresponding uncertainties. 

For each cluster, the field-star contamination was re¬ 
moved from the cluster GMD, and Monte Garlo realiza¬ 
tions of this background-subtracted cluster GMD were 
created. Monte Carlo isochrone realizations, based on 


Fig. 14.— Locations of the 30 fit clusters in the outer Galaxy, 
with the Sun at the origin. The clusters are shown as filled circles, 
and the color denotes the degree of reddening, as indexed by the 
legend in the bottom left corner. The sizes of the circles correspond 
to the number of cluster CMD members (Table [^. as indicated by 
the legend in the upper right corne r. The approximate location of 
the Perseus Arm lIHeid et al.ll2Ql^ is shown as the grey stripe. 

the PARSEC isochrones, were created at different ages, 
distance moduli and color excesses to compare to each 
cluster CMD. The isochrone parameters which yielded 
the lowest color-magnitude distance-based mean were 
adopted as the cluster parameters. Of the original sam¬ 
ple of 31 clusters, 24 were fit for all three parameters. 
The remaining seven were either faint or sparse. By fix¬ 
ing age, distance and reddening estimates were fit for six 
of the seven. The mean uncertainties of the 30 clusters 
of log(age), distance, and E{B — V) were 0.13, 340 pc, 
and 0.14 mag, respectively. 

To test the accuracy of the fitting technique, synthetic 
clusters were created and fit. These clusters were gen¬ 
erated at 12 ages, ranging in log(age) from 7.3 to 9.6, 
for 50, 100, and 200 numbers of stellar members. For 
nearly all of the synthetic clusters, the input parame¬ 
ters were recovered to within their la uncertainties. The 
clusters with 50 members were the most difficult to fit 
and showed the largest scatter from their input values, 
though a strong majority were still fit to within Icr of 
their input parameters. 

The derived properties of the 30 clusters that were fit 
were compared to values found in recent published stud¬ 
ies and revealed no biases or trends. 

This sample of clusters spans wide ranges in longitude 
and distance. It also probes the Perseus Spiral arm, as 
well as its foreground and background inter arm regions. 
The cluster distribution is well suited to probe the prop¬ 
erties of the large-scale Galactic magnetic field in the 
outer Galaxy and to test whether the field is affected by 
the presence of a spiral arm. 
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TABLE 1 

Open Cluster Sample and Stellar Radial Density Profile Parameters 


Name 

1 

(deg) 

b 

(deg) 

Number of 
Stars 

(arcsec) 

Background 
(Stars arcmin”^) 

Metallicity 

(Z) 

Metallicity 

Ref.l’ 

Berkeley 60 

118.846 

-1.643 

130 ±50 

180 ±29 

6.21 ±0.07 

0.019 

— 

King 1 

119.763 

1.693 

220 ±50 

161 ±15 

6.96 ±0.07 

0.019 

— 

NGC 559 

127.194 

0.746 

360 ±80 

275 ±25 

6.32 ±0.07 

0.019 

— 

NGC 663 

129.467 

-0.941 

520 ±100 

318 ±25 

6.72 ±0.08 

0.019 

— 

NGC 869 

134.633 

-3.752 

1020 ±150 

455 ±27 

5.47 ±0.09 

0.019 

— 

King 5 

143.773 

-4.276 

220 ±50 

177 ±15 

4.07 ±0.06 

0.008 

(Friel 19951 

NGC 1245 

146.654 

-8.908 

430 ±70 

280 ±18 

3.27 ±0.05 

0.017 

f Jacobson et al. 2011) 

King 7 

149.799 

-1.021 

490 ±70 

201 ±11 

5.02 ±0.06 

0.019 

— 

Berkeley 12 

161.677 

-1.992 

170 ±40 

145 ±14 

4.64 ±0.06 

0.019 

— 

Berkeley 14 

162.873 

0.707 

170 ±50 

174 ±19 

5.21 ±0.06 

0.019 

■— 

NGC 2126 

163.236 

13.140 

54 ±24 

124 ±21 

2.11 ±0.04 

0.019 

— 

Berkeley 18 

163.648 

5.050 

410 ±100 

410 ±40 

3.91 ±0.07 

0.019 

— 

Berkeley 70 

166.902 

3.585 

200 ±70 

300 ±40 

4.11 ±0.06 

0.019 

fCarrera 20121 

NGC 1857== 

168.444 

1.222 

360 ±100 

390 ±40 

4.97 ±0.07 

0.005 

— 

NGC 2099 

177.636 

3.092 

1030 ±140 

448 ±24 

5.12 ±0.08 

0.019 

I'Pancino et al. 2010) 

NGC 2158 

186.634 

1.781 

890 ±80 

192 ±6 

5.87 ±0.07 

0.010 

f Jacobson et al. 2011) 

Basel 11b 

187.442 

-1.114 

32 ±21 

90 ±24 

4.88 ±0.06 

0.019 

— 

NGC 2266 

187.778 

10.304 

88 ±25 

94 ±9 

2.53 ±0.04 

0.007 

(Reddv et al. 2013) 

NGC 2141 

198.044 

-5.810 

440 ±60 

218 ±12 

3.59 ±0.05 

0.019 

(Jacobson et al. 2009^ 

NGC 2420 

198.107 

19.634 

360 ±50 

225 ±12 

1.45 ±0.04 

0.012 

(Jacobson et al. 2011) 

Trumpler 5 

202.807 

1.018 

1510 ±140 

434 ±16 

4.59 ±0.08 

0.008 

(Garrera et al. 200?') 

NGC 2355 

203.390 

11.803 

230 ±50 

246 ±21 

2.21 ±0.04 

0.016 

(Jacobson et al. 2011) 

NGC 2112 

205.872 

-12.615 

420 ±100 

470 ±50 

2.20 ±0.06 

0.015 

IBrown et al. 1996^ 

Berkeley 32 

207.952 

4.404 

280 ±70 

283 ±28 

3.82 ±0.06 

0.0095 

(Garrera V Pancino 2011) 

Collinder 110 

209.649 

-1.927 

530 ±120 

460 ±40 

4.66 ±0.08 

0.019 

(Garrera et al. 2007) 

NGC 2262 

210.573 

-2.099 

190 ±50 

151 ±14 

4.42 ±0.06 

0.019 

— 

NGC 2324 

213.447 

3.297 

210 ±50 

198 ±19 

4.39 ±0.06 

0.0128 

fBraeafflia et al. 2008) 

M 67 

215.696 

31.896 

450 ±90 

520 ±40 

1.03 ±0.04 

0.019 

(Jacobson et al. 2011) 

Berkeley 39 

223.461 

10.095 

210 ±40 

186 ±15 

2.34 ±0.04 

0.014 

(Garrera et al. 200?') 

Haffner 10 

230.799 

1.011 

250 ±60 

185 ±16 

5.49 ±0.07 

0.019 

— 

NGC 2425 

231.501 

3.289 

150 ±50 

184 ±23 

4.28 ±0.06 

0.0135 

fJacobson et al. 2011) 


^ Pi-SGLU signifies the radius of the inner region, which is equal to 1.5 times the Gaussian width of the cluster radial density 
profile. 

^ For clusters where no spectroscopic metallicity measurement was found, no reference is listed. 

^ For cluster NGC 1857, a metallicity of Z = 0.005 was assumed instead of solar metallicity. 
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TABLE 2 

Derived Cluster Properties 


Name 

Number of ^ 

Log (age) 


E(J-if) 

Distance 

E(B - VY 


CMD Stars 

[log(yrs)] 

(mag) 

(mag) 

(pc) 

(mag) 


Berkeley 60 

56 ±16 

7.88 ±0.13 

14.46 ±0.22 

0.59 ±0.07 

6160 1680 

1.11 

King 1 

135 ±17 

9.22 ±0.15 

11.74 ±0.21 

0.42 ±0.07 

1880 l?oo 

0.80 

NGC 559 

218 ±26 

9.20 ±0.13 

11.80 ±0.30 

0.23 ±0.07 

2070 1840 

0.42 

NGC 663 

260 ±30 

7.25 ±0.13 

12.96 ±0.27 

0.42 ±0.07 

3300 1^0 

0.80 

NGC 869 

580 ±40 

7.25 ±0.13 

12.30 ±0.30 

0.26 ±0.09 

2560 l«o 

0.50 

King 5 

114 ±15 

9.15 ±0.12 

11.69 ±0.20 

0.38 ±0.08 

1870 IJ^^J 

0.71 

NGC 1245 

243 ±21 

9.18 ±0.12 

12.39 ±0.20 

0.09 ±0.07 

2900 l^?g 

0.17 

King 7 

223 ±20 

8.03 ±0.12 

13.59 ±0.22 

0.78 ±0.07 

3820 l«o 

1.46 

Berkeley 12 

57 ±12 

9.53 ±0.11 

12.99 ±0.18 

0.35 ±0.08 

3440 1820 

0.66 

Berkeley 14 

50 ±14 

O.SO*" 

14.00 ±0.19 

0.25 ±0.07 

5700 1880 

0.47 

NGC 2126 

32 ±8 

9.45'" 

10.10 ±0.40 

0.09 ±0.08 

1010 l?io 

0.17 

Berkeley 18 

177 ±27 

9.63'" 

13.78 ±0.21 

0.26 ±0.08 


0.50 

NGC 1857 

166 ±28 

8.50 ±0.12 

12.24 ±0.21 

0.31 ±0.07 

2470 

_240 

0.59 

NGC 2099 

640 ±40 

8.78 ±0.11 

10.75 ±0.17 

0.17 ±0.07 

1320 1^1° 

0.33 

NGC 2158 

372 ±23 

9.25 ±0.12 

13.23 ±0.18 

0.25 ±0.07 

3990 

oyyu _33o 

0.47 

Basel lib 

29 ±8 

9.00'" 

9.72 ±0.29 

0.21 ±0.07 

Q-| (3 +120 

oJ-0 _;i^QQ 

0.40 

NGC 2266 

44 ±8 

9.15 ±0.12 

12.40 ±0.20 

0.14 ±0.08 

2860 IjOO 

0.26 

NGC 2141 

181 ±18 

9.25 ±0.14 

12.99 ±0.29 

0.29 ±0.07 

3520 1820 

0.54 

NGC 2420 

147 ±14 

9.22 ±0.16 

12.10 ±0.25 

0.07 ±0.08 

2550 1880 

0.14 

Trumpler 5 

620 ±40 

9.53 ±0.12 

12.80 ±0.18 

0.35 ±0.07 

3150 1^00 

0.66 

NGC 2355 

115 ±15 

9.15 ±0.14 

11.43 ±0.23 

0.10 ±0.07 

1850 l?io 

0.19 

NGC 2112 

273 ±27 

9.53 ±0.13 

9.81 ±0.22 

0.30 ±0.07 

810 10° 

0.57 

Berkeley 32 

108 ±20 

9.45'" 

13.10 ±0.30 

0.01 ±0.08 

4150 16^0 

0.02 

Collinder 110 

320 ±40 

9.20 ±0.16 

11.91 ±0.24 

0.20 ±0.08 

2230 l^^O 

0.38 

NGC 2262 

92 ±14 

9.07 ±0.15 

12.54 ±0.21 

0.33 ±0.08 

2820 1810 

0.61 

NGC 2324 

105 ±16 

8.70 ±0.14 

13.30 ±0.21 

0.14 ±0.07 

4320 1«0 

0.26 

M 67 

327 ±24 

9.47 ±0.15 

9.15 ±0.23 

0.04 ±0.07 

670 l^p 

0.07 

Berkeley 39 

70 ±12 

9.50'" 

13.21 ±0.20 

0.11 ±0.08 

4200 l«o 

0.21 

Haffner 10 

99 ±17 

9.50 ±0.13 

13.13 ±0.18 

0.23 ±0.07 

3850 I860 

0.42 

NGC 2425 

85 ±16 

9.18 ±0.14 

12.88 ±0.26 

0.25 ±0.07 

3400 1^« 

0.47 


^ The numbers of CMD stars are the total number of star counts of the cluster CM P overdensity 
images, and represent the counts the Poisson images are drawn from. See Section[3]4|for discussion. 
^ Uncertainties in log(age) are not listed, as these clusters w ere only fit for distance m odulus and 
color excess, and their ages were fixed at the values listed in IKharchenko et al.l 1I2013IV 
^ The uncertainties of E(S — V) are equal to the uncertainties of E(J — iC)/0.53, following the 
relation between E(B — V) and E(J — K). 
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TABLE 3 

Cluster Parameters Found in Recent Studies 


Name 

Log (age) 
[log(yrs)] 

D 

(pc) 

E(B-V) 

(mag) 

Agreement' 

Berkeley 60 

8.25 

4468 

0.92 

ccb 

Berkeley 60 

8.2 

4365 

0.86 

ccb 

Berkeley 60 

8.4 

3299 

0.958 

ddb 

Berkeley 60 

8.2 

2089 

0.37 

cdd 

King 1 

9.45 

2060 

0.625 

bab 

King 1 

9.6 

1080 

0.76 

cda 

King 1 

9.2 

1900 

0.7 

aaa 

King 1 

9.7 

1659 

0.625 

dbb 

NGC 559 

8.35 

2430 

0.82 

dbd 

NGC 559 

8.8 

2170 

0.68 

dab 

NGC 559 

8.6 

2291 

0.81 

dac 

NGC 559 

7.7 

6309 

0.62 

ddb 

NGC 559 

— 

1200 

— 

-c- 

NGC 559 

9.08 

1300 

— 

ac- 

NGC 559 

8.8 

2200 

0.6 

dab 

NGC 663 

7.65 

2520 

0.63 

dbb 

NGC 663 

7.3-7.4b 

2420 

0.8 

aba 

NGC 663 

7.4 

2089 

— 

be- 

NGC 663 

7.3 

2469 

0.75 

aba 

NGC 663 

7.3-7.4 

2100 

0.83 

aca 

NGC 663 

7.86 

1718 

— 

dd- 

NGC 663 

7.08-7.4 

2818 

0.8 

aba 

NGC 663 

7.5 

2100 

0.7 

bca 

NGC 869 

7.14 

2290 

0.55 

aaa 

NGC 869 

7.11 

2269 

0.54 

baa 

NGC 869 

— 

— 

0.52 

—a 

NGC 869 

7.11 

2344 

0.56 

baa 

NGC 869 

7.1 

2188 

— 

ba- 

NGC 869 

7.1 

2904 

0.58 

baa 

NGC 869 

6.91 

2025 

— 

cb- 

NGC 869 

7.28 

2300 

0.521 

aaa 

King 5 

9.1 

2230 

0.67 

aba 

King 5 

9 

1900 

0.82 

baa 

King 5 

9 

1905 

0.94 

bab 

King 5 

8.9 

— 

— 

c- 

King 5 

9.09 

2200 

0.67 

aba 

NGC 1245 

9.03 

2818 

0.24 

baa 

NGC 1245 

9.04 

3010 

0.05 

baa 

NGC 1245 

9.02 

2800 

— 

ba- 

NGC 1245 

8.95 

3019 

0.29 

baa 

NGC 1245 

9.03 

— 

— 

b- 

NGC 1245 

9.16 

2211 

0.27 

aca 

NGC 1245 

9.025 

3000 

0.25 

baa 

King 7 

8.7 

— 

— 

d- 

King 7 

8.78-8.9 

2200 

1.25 

ddb 

King 7 

8.6 

2440 

1.25 

ddb 

King 7 

8.85 

2628 

1.249 

deb 

Berkeley 12 

9.2 

3801 

0.8 

dba 

Berkeley 12 

9.6 

3162 

0.7 

aaa 

Berkeley 12 

9.6 

3300 

0.7 

aaa 

NGC 1857 

9 

1400 

0.13 

ddd 

NGC 1857 

8.0-8.25 

5750 

0.38-0.6 

dda 

NGC 1857 

8.2 

1545 

0.97 

ede 

NGC 1857 

8.67 

3299 

0.5 

bda 

NGC 2099 

8.34-8.51 

— 

— 

d- 

NGC 2099 

8.69 

1490 

0.227 

aba 

NGC 2099 

— 

1905 

0.21 

-da 

NGC 2099 

00 

bo 

1995 

0.23 

ad a 

NGC 2099 

8.6-8.72 

1148-1202 

0.36 

bba 

NGC 2099 

8.55 

1400 

0.35 

caa 

NGC 2158 

9.28 

3944 

0.42 

aaa 

NGC 2158 

9.28 

— 

— 

a- 

NGC 2158 

9.3 

3600 

0.55 

aba 

NGC 2158 

— 

4068 

0.43 

-aa 

NGC 2158 

9.2 

5012 

0.4 

aca 

NGC 2158 

9.33 

4770 

0.333 

acb 

NGC 2266 

9 

2855 

0.21 

baa 

NGC 2266 

9.08 

2800 

0.17 

aaa 

NGC 2266 

8.94 

— 

— 

b- 

NGC 2266 

8.9 

3758 

0.1 

cdb 

NGC 2266 

9.265 

3311 

0 

abb 

NGC 2141 

9.1-9.28 

4090-4370 

0.36-0.45 

abb 

NGC 2141 

9.4 

— 

— 

b- 


Ref. 


B_ukowieck^^aL ^01^ 


JOiarchenk^^^aL ^0^) 

Hasegawa et al. (2008) 
^acieiewsk^^^JiedzielsM ^00^ 

JOiarchenk^^^aL 

^aciej^^^^^^^^dzielsM ^00^ 

Jennens &: Heifer f l975) 
Grubissich (1975) 
^^indo^ Tl96^) 
^harchenk^^^aL 
^ukowieck^^aL ^OH) 
^ande^^^al^ ^00^ 
Fabregat &: Capill^ ^OOS) 
( 2001 ) 

PheIp^^J^e^ ^^94‘l 

^harchenk^^^al^ 

^a^ne^^^Ja^lOT (gSo^ 
^Bragg_&_Ken^on (^200^ 

Capin£^^Fabr^at^ ^02 ) 

^al^she^ (T9^) 

( 2001 ) 

^arraro_^^allen^ (2000) 
JCharchenk^^^T^^^^) 

^urke_et_^ (2004) 
^ubramaniam^O^) 
^alari^^ak (2004 ) 
^^^dros^ ^OOl) 
JCharchenk^^^^L (^^^) 

^urga£a^^^^^97) 

J^harchenk^^^al^ ^^^) 
jlasegaw^^^al^ QOO^ 

J^harchenk^^^^^^^^) 
^asowsk^^al^ ^0^) 
^ujath^^^al^ 

^^Tadros ^^O^^ "^ 
JCharchenk^^^L (^^^) 

JJartmai^^ak^^OOS) 

JCharchenk^^^^A^^S) 

^alari^^al^ 
^arrar^^^^L ^0^) 
^rocholsk^^^amjedim ^002) 
^adro^ (2001) 
JCharchenk^^^a^ (^^^) 
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Finding Open Cluster Properties 


TABLE 3 — Continued 


Name 

Log (age) 
[log(yrs)] 

D 

(pc) 

E(B-V) 

(mag) 

Agreement' 

NGC 2141 

9.4 

3800 

0.4 

bab 

NGC 2141 

9.4 

4200 

0.35 

bbb 

NGC 2141 

9.245 

4364 

0.312 

abb 

NGC 2420 

9.3 

2480 

0.04 

aaa 

NGC 2420 

9.3 

2542 

0.04 

aaa 

NGC 2420 

9.34 

2443 

0.05 

aaa 

NGC 2420 

9.3 

2449 

0.05 

aaa 

NGC 2420 

9.365 

2880 

0.01 

aaa 

Trumpler 5 

9.7 

2400 

0.5 

bcb 

Trumpler 5 

9.54-9.6 

2818-3076 

0.6-0.7 

aaa 

Trumpler 5 

9.45 

3100 

0.64 

aaa 

Trumpler 5 

9.75 

— 

— 

b- 

Trumpler 5 

9.61 

3019 

0.58 

aaa 

Itumpler 5 

9.1 

2958 

0.58 

daa 

Trumpler 5 

9.5 

2753 

0.625 

aba 

Trumpler 5 

9.5 

2900 

0.6 

aaa 

NGC 2355 

8.9 

1985 

0.3 

baa 

NGC 2355 

8.9 

— 

— 

b- 

NGC 2355 

9 

1650 

0.16 

baa 

NGC 2355 

8.98 

1915 

0.112 

baa 

NGC 2355 

8.9 

2128 

0.187 

bba 

NGC 2112 

9.23 

940 

0.6 

cba 

NGC 2112 

9.45 

813 

0.6 

aaa 

NGC 2112 

9.6 

750 

0.6 

aaa 

NGC 2112 

9.315 

977 

0.625 

bba 

Collinder 110 

9.08-9.23 

1949-2187 

0.52-0.58 

aab 

Collinder 110 

9.15 

1950 

0.5 

aba 

Collinder 110 

9.22 

2362 

0.416 

aaa 

NGC 2262 

9 

3600 

0.55 

aca 

NGC 2262 

8.995 

2511 

0.625 

aaa 

NGC 2324 

8.65 

3800 

0.25 

aba 

NGC 2324 

8.83 

— 

— 

a- 

NGC 2324 

8.8 

4169 

0.17 

aaa 

NGC 2324 

8.68 

3842 

0.239 

aba 

NGC 2324 

8.8 

4400 

0.1 

aab 

M 67 

9.05 

722 

0.24 

cab 

M 67 

9.6 

823 

0.04 

aba 

M 67 

9.56-9.66 

795 

0.038 

aba 

M 67 

— 

766 

0.038 

-ba 

M 67 

9.6 

832 

0.04 

aca 

M 67 

9.51 

870 

0 

aca 

M 67 

9.535 

890 

0.05 

aca 

Haffner 10 

9.2-9.4 

3700 

0.55 

baa 

Haffner 10 

9.4 

3100-4300 

0.41-0.64 

aaa 

Haffner 10 

9.305 

4873 

0.5 

bca 

NGC 2425 

9.56 

3357 

0.175 

cac 

NGC 2425 

9.34 

3550 

0.21 

bab 

NGC 2425 

9.4 

2900-3800 

0.29 

bab 

NGC 2425 

9.34 

4330 

0.21 

bcb 


Ref. 


Carraro et al. (^01) 

^Rosvi^ 719951 
I^harchenk^^e^^^ (2013) 
Sharma et al. (2006 ^ 
Mermilliod &: Mayoi^ 72007 ) 
Salaris et al. (2004) 
^rocholsk^^^Saraied^ (2003 
^harche^o et aL (^13 ) 
^Piatti^^^^ f 200i) 
_Don^l^^aL 

(2009‘l 

Salaris et al. (20^4) 
KaluznvTi998) 
^Tadro^ (2 001) 
Kharchenko et al. (2013 ) 
'^~ "Perren et al. (201^ " 
^Dias_et_^L (201? ) 
Salaris et al. (^04) 
^oubirai^^^_^^ ) 

Tadross ( '^00j3 
Kharchenko et al. (2013 ) 
Carraro et al. . 

Tadross (20^5 
Richtler &; Kaluzny (^989) 
Kharchenko et al. (2013 ) 
^ragaglia &: T<^i (2003 ) 
^awsoi^_^J^na_ (1998 ) 
Kharchenko et al. (2013 ) 
^ ^arraro e" a !l. r2005 ) 
Kharchenko et ahT^ 13) 
Piatti et al. 72004b) 
Salaris et al. (2004) 
Kyeong et al. (2001) 
JOiarchen^jst al. ( 2013 ) 
Perren et al. (20l5) 

Dias et a l. (2012) 
Balagner-Nune z et ah (20pJ) 
VandenBerg Stetson f%6 4) 
Laugalys et al. (2004^^ 
Sandquis t (2004) 
B onatto Sz B ica (2003) 
K harchenko et al. ( 2013 ) 
Vaz^uez^^ al. 
Pie trukowicz et al, 
Kh archenko et ah 
Ilaseg^v^^^^ 

Moitinh^et^ 
^intrukowic^e^h 

Kharchenko et al. 



^ Quality of agreement within 1, 2, 3, or > Str: a, b, c, d, respectively. 
^ For properties reported as a range, the average of the range was 
used for comparison. 






















































